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of mixed Gaussian states. A comparison between mutual information and entanglement of formation for symmetric states is 
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1 Introduction 

Quantum information with continuous variable systems is 
rapidly developing and appears to yield very promising 
perspectives concerning both experimental realizations and 
general theoretical insights. In such a context, Gaussian 
states play a prominent role, both in view of their concep- 
tual importance and for their relevance in experimental ap- 
plications, and have attracted most of the attention of the 
researchers in the field 1 1 1. They are the easiest states that 
can be created and controlled in the laboratory |2 1, and have 
been successfully exploited in quantum cryptography 1 3 1 
and quantum teleportation protocols |4|. Moreover, they 
are possible candidates for continuous variable quantum 
computation processing |5 1 1. 

In particular, two-mode Gaussian states aroused great 
interest in later years, being the simplest prototype of a 
continuous variable bipartite system. As for the theory, 
the qualitative characterization of the entanglement of two- 
mode Gaussian states has been fully developed by deter- 
mining the necessary and sufficient criteria for their sepa- 
rability EE). A computable quantitative characterization 
of entanglement for such states is available as well, being 
provided by negativity 1 8 1 and, in the symmetric instance, 
by the entanglement of formation I9l ll0l . 

Due to the interaction with the environment, any pure 
quantum state involved in some quantum information pro- 
cess evolves into a mixed state. Therefore, another property 
of crucial interest in Quantum Information Theory is quan- 
tifying the degree of mixedness of a quantum state. Let 
us briefly recall that the degree of mixedness of a quantum 



state g can be characterized either by the Von Neumann en- 
tropy Sv{q) or by the linear entropy Sl{q)- Such quanti- 
ties are defined as follows for continuous variable systems: 

Sv{q) = -Trig lug), (1) 

SUg) = 1- ^Tig^) = l-^l{g), (2) 

where /i = Tr (g^) denotes the purity of the state g. The 
linear entropy of an arbitrary n-mode Gaussian state can be 
easily computed, whereas the evaluation of the Von Neu- 
mann entropy requires, in general, a more involved tech- 
nical procedure. Making use of the Von Neumann en- 
tropy is however preferable, as it allows for a finer and 
more precise characterization of mixedness and correla- 
tions for multi-mode Gaussian states. In fact, the Von 
Neumann entropy is additive on tensor product states, un- 
like the linear entropy. Moreover, we will show that for 
two-mode Gaussian states the Von Neumann entropy de- 
pends on two symplectic invariants, while the linear en- 
tropy in completely determined by only one invariant. Fi- 
nally, the knowledge of the Von Neumann entropy of a 
generic two-mode state allows to obtain the mutual infor- 
mation I{g) = Sv{gi) + Sv{g2) — Sv{q) (here g^ is the 
reduced density matrix of subsystem i), which quantifies 
the total amount of correlations (quantum plus classical) 
contained in a state 1 1 1 1 . 

The Von Neumann entropy of a n-mode Gaussian state 
has been determined in some remarkable works on the ca- 
pacity of Gaussian channels by Holevo et al. 1121 . In this 
letter we elucidate the mathematical structure introduced 
by these Authors, recasting it in a form that is conceptu- 
ally simple and physically useful for applications and for 
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explicit experimental schemes. To do so, we provide a 
thorough and detailed analysis of the basics of the sym- 
plectic framework restricting to the two-mode case, pre- 
senting a simple and transparent procedure to evaluate the 
Von Neumann entropy and the mutual information of two- 
mode Gaussian states. Our approach is based on the ex- 
plicit determination of the global symplectic invariants. 
The methodology we introduce naturally highlights the role 
played by the symplectic eigenvalues in characterizing the 
correlations encoded in multipartite systems. 

2 Two-mode Gaussian states 

Let us consider a two-mode continuous variable system, 
described by the Hilbert space H = 7^1(8)7^2 resulting from 
the tensor product of the Fock spaces Ti^'s. We will call ak 
the annihilation operator acting on the space TCk- Likewise, 



[ak 



:)/^/2 and 



Pk 



,)/\/2 are the 



quadrature phase operators of the mode k, the correspond- 
ing phase space variables being Xk and pk- The set of Gaus- 
sian states is, by definition, the set of states with Gaussian 
characteristic functions and quasi-probability distributions. 
Therefore, a Gaussian state is completely characterized by 
its first and second statistical moments, that is, respectively, 
by the vector of mean values X = ((xi), (pi), {£2), {P2)) 
and by the covariance matrix cr 



(3) 



First moments will be unimportant to our aims, and we will 
set them to zero (as it is always possible by means of a 
local unitary transformation) without any loss of generality 
for our results. For simplicity, in what follows cr will refer 
both to the Gaussian state and to its covariance matrix. It is 
convenient to express cr in terms of the three 2x2 matrices 
a, (3, 7 

'a 7 



7^ /3 



Let us define two further submatrices of cr 
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The privileged role played by 6 and e in characterizing the 
action of global symplectic operations on cr will become 
clear in the following. 

Positivity of g and the commutation relations for quadra- 
ture phase operators impose the following constraint ensur- 
ing that cr be a bona fide covariance matrix li6j| 



-O > , 



(6) 



where ft is the usual symplectic form 





















(7) 



Inequality (|6|l is a useful way to express the Heisenberg 
uncertainty principle. 

In the following, we will make use of the Wigner quasi- 
probability representation W{xi,pi), defined as the Fourier 
transform of the symmetrically ordered characteristic func- 
tion |13|. In the Wigner phase space picture, the tensor 
product H = Hi <8i 7^2 of the Hilbert spaces Hi's of the 
two modes results in the direct sum F = Fi r2 of the re- 
lated phase spaces F^'s. As a consequence of the Stone- Von 
Neumann theorem, a symplectic transformation acting on 
the global phase space F corresponds to a unitary operator 
acting on H \ 14|. In what follows we will refer to a trans- 
formation Si — Si (B S2, with each Si G S'p(2,R) acting on 
Fi, as to a "local symplectic operation". The corresponding 
unitary transformation is the "local unitary transformation" 
Ui = Ui® U2, with each Ui acting on Hi- Inequality ^ is 
then a constraint on the S'p(2,m) © <5'P(2,r) invariants (6| 

Det a + Det ;3 + 2 Det 7 < i + 4 Det cr . 

The Wigner function of a Gaussian state, written in terms 
of the phase space quadrature variables, reads 



W{X) 



7r\/E)et cr 



(8) 



where X stands for the vector {xi,pi,X2,P2) G T. In 
general, the Wigner function transforms as a scalar under 
symplectic operations, while the covariance matrix cr trans- 
forms according to 

cr S'^a-S , S e 5p(4^R) . 

As it is well known, for any covariance matrix cr there ex- 
ists a local canonical operation Si = SiG) S2 which brings 
cr in the "standard form" (Ts/ Q 



Si crSi — 



a 





Cl 







a 





C2 


Cl 





h 








C2 





b 1 



(9) 



where a, 5, ci, C2 are determined by the four local symplec- 
tic invariants Deter = [ab — c\){ah — C2), Detct = a?. 
Dot /3 — b^, and Det 7 = C1C2. Therefore, the coefficients 
of the standard form corresponding to any covariance ma- 
trix are unique (up to a common sign flip of the c^'s). 
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3 Determination of the Von Neumann 
entropy 

To proceed, let us first note that the purity /i (and therefore 
the linear entropy Sl) of a Gaussian state can be easily 
computed. In fact, the trace of a product of operators cor- 
responds to the integral of the product of their Wigner rep- 
resentations (when existing) over the whole phase space. 
Using the Wigner representation W of g and taking into ac- 
count the proper normahzation factors, for a n-mode Gaus- 
sian state we get 



1 



2'V Deter 



(10) 



Eq. ilO\ implies that a Gaussian state cr is pure if and only 
if Deter = 1/22". 

For single-mode systems, the Von Neumann entropy can 
be easily computed as well. Neglecting first moments, any 
single-mode Gaussian state g can in fact be written as 

g = Ss,n{r, ip)unSl^{r, ip) , (11) 

where Un is a thermal state of mean photon number n 



fe=0 



\k){k\ 



(12) 



e^^'^aT^) is the 



single-mode squeezing operator. Being unitary, the lat- 
ter does not affect the values of the traces in Eqns. Q- 
Q, computed on the diagonal density matrix Un given by 
Eq. M2\ . One has then 



Sv{q) = n\n 
2^ 



1 



1 



2n + 1 2V Det cr 



(13) 



hi 



ln(rl + 1) 



-In 



1 



(14) 



Eq. ( I14> . first achieved in Ref. 1 15 1, shows that for single- 
mode Gaussian states the Von Neumann entropy is a mono- 
tonically increasing function of the linear entropy, so that 
Sv and Sl yield the same characterization of mixedness. 
In fact, both Sy and Sl are fully determined by the same 
symplectic invariant Det er. As we will now see, this is no 
longer true for two-mode Gaussian states. 

To find an expression for the Von Neumann entropy of 
a generic Gaussian state of a two-mode system, we must 
find a general expression for the state analogous to that pro- 
vided by Eq. (II 1> for a single-mode system. Neglecting 



first moments, this amounts to determine the most general 
parametrization of the covariance matrix, which is provided 
by the following lemma. 

Lemma 1 An arbitrary two-mode covariance matrix cr can 
be written as 

cT^A^u^^A, (15) 

where i/„^ = i>'„__i/2 ® is the covariance ma- 

trix of a tensor product of single-mode thermal states with 
average photon number flip = — 1/2 in the two modes 



while 



iv„ = diag(n_,7i_,7i+,7i+) , (16) 



A = Sioc{ri,r2)R{0StUr)R{v)Si (17) 



is a symplectic operation belonging to 5p(4_R). Transfor- 
mation A is made up by a local operation Si, two rotations 
with 
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cos <l 



I 



(18) 



a two-mode squeezing Stmij) = diag(e'', e"'', e"'', e'') 
and a local squeezing 5/oc(fi,f2) — 5'sm (ri,0) ® 
Ssm (''2 5 0), resulting from the direct product of two single- 
mode squeezing operators with null phase. Note that 
Stjn{r) = Sioc{r, —r), so that the only global (nonlocal) 
operations in the decomposition of Eq. ( I17> are the two ro- 
tations. We note that an equivalent decomposition has been 
recently demonstrated for generic multimode pure Gaus- 
sian states 1 16 1: the authors have shown that for any decom- 
position of a multimode pure Gaussian state with respect to 
a bipartite division of the modes, the state can always be 
expressed as a product state involving entangled two-mode 
squeezed states and single-mode local states at each side. 

Proof. In order to prove the statement expressed 
by Eq. il5\ . we consider the equivalent expression 
A^^'^crA^^ = and show that it is realized by some 
A, n_, n+ for any given cr. 

First, we choose Si to bring cr to its standard form, given 
by Eq. (|9j. We then apply to a-gf the rotated two-mode 
squeezing R{T]y^St„i{ry^, 
to the form 






V 




m 


c 



taking the covariance matrix 



which is convenient, due to the invariance of the submatrix 
6 — diag(s, s), see Eqns. (|5j, under two-mode rotations 
of the form Eq. ( I18> . The second rotation leaves 
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S unchanged and can be chosen to make c null, yielding a 
state of the form 

/ s ^ 
m' 
s 
\ On'/ 

which can be finally put in the desired form i>'„__i/2 © 
i'n+-i/2 by means of the local squeezing Sioc{ri, ^2). □ 

Lemma 1 introduces an equivalence relation on the set 
of Gaussian states, associating to any Gaussian state cr a 
product of thermal states I'n^, by means of the correspon- 
dence defined by Eq. (I15> . The quantities rizp are known 
as the symplectic eigenvalues of cr, while transformation A 
performs a symplectic diagonalization | IT'S* 171. We note 
that the decomposition of the symplectic operation A given 
by Eq. (tTTt is the particular two-mode case of the general 
decomposition of a symplectic operation 1181 . 

Let us now focus on the quantity 

A((t) = Dcta + Det/3 + 2Det7, (19) 

where a, f3, and 7 are defined as in Eq. (0}. We have: 

Lemma 2. A((t) is invariant under the action of the sym- 
plectic transformation A defined by Eq. (I15> . 

Proof. A(cr) is clearly invariant under local operations, 
such as Si, Sloe and Stm- As for the non local rotations 
which enter in the definition of A, let us notice that they act 
on covariance matrices of the following form 



The solution of the system yields 



u 

j 

V 





V 



k 



J 

w 




k 


^ J 



for which one has 



A(c^) = Tr (<5e) 



Such an expression is manifestly invariant under the action 
of identical rotations R{4>) on the submatrices 6 and 7, see 
Eqns. (|5} and ([Ts}. □ 

The quantity Det cr is obviously invariant as well under 
the action of A since, for any symplectic transformation S, 
one has Det S* = 1. Exploiting the invariance of Deter 
and A(cr) one can determine the symplectic eigenvalues 
which characterize a generic Gaussian state cr, accord- 
ing to Eq. (ITsl 



Deter — Detv>n^ = n^n^ , 



A(cr) = A{VnJ ^nl+nl 



(20) 



nzp(cr) 



A(cr) =F v/A(cr)2 -4Dctcr 



Note that Ineq. (|6j is equivalent to 



> — , 



(21) 



(22) 



whereas the necessary and sufficient criterion for a state to 
be pure reads n_ = n+ = 1/2 (one can easily show that it 
is equivalent to Det cr = 1/16). 

Knowledge of the symplectic eigenvalues and of the as- 
sociated mean thermal photon numbers allows finally to de- 
termine the Von Neumann entropy S'y(cr) of an arbitrary 
two-mode Gaussian state cr. We have: 
Proposition 1. The Von Neumann entropy 5v(cr) of an ar- 
bitrary two-mode Gaussian state cr equals the one of the 
tensor product of thermal states associated to or via 
the correspondence established by Eq. il5\ . and its expres- 
sion reads 



Sv{a)^f[n-{a)] + f[n+{a)] 



(23) 



l)Hx-^). (24) 



with rizp given by Eqns. ( 12 H and 

f{x) = {x + ^)\n{x + ^)-{x 

Proof The symplectic operation A described by Eq. J17> 
corresponds to a unitary transformation in the Hilbert space 
H which cannot affect the value of the trace appearing in 
the definition of Sy, according to Eq. Q. Therefore, ex- 
ploiting Eq. il4i and the additivity of the Von Neumann 
entropy for tensor product states, one obtains Eq. ( I23t . □ 

We have shown that the Von Neumann entropy of a two- 
mode Gaussian state cr depends on the two invariants A(cr) 
and Det cr, whereas the purity of cr is completely deter- 
mined by Det cr alone, just as in the single-mode case. 
This implies that the hierarchy of mixedness established by 
the Von Neumann entropy on the set of Gaussian states dif- 
fers, in the two-mode case, from that induced by the lin- 
ear entropy. States may exist with a given linear entropy, 
i.e. with a given Deter, but with different Von Neumann 
entropies, i.e. with different A(cr)'s. The Von Neumann 
entropy thus provides a richer characterization of the state's 
lack of information. 

4 Symplectic eigenvalues, mutual infor- 
mation and correlations 

The mutual information /(cr) of a Gaussian state cr is de- 
fined as 

I{(t)^Sv{(ti) + Sv{(T2)-Sv{(t) , (25) 
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where ai stands for the reduced single-mode state obtained 
by tracing over subsystem j ^ i. Knowledge of S'y(cr) 
leads to the following: 

Proposition 2. The mutual information of an arbi- 
trary two-mode Gaussian state is 

= f{a) + f{b) - /[n_(<T)] - f[n+{cT)] , (26) 

where a — %/ Det a. , b — \J Det /3 , andf(x) is the same 
as in Eq. A241 . 

Proof. Let us consider the reduction of cr to its standard 
form (Tsf, defined by Eq. (|9j. The matrix elements a and 
6 of (Tsf are easily recovered from a generic cr, because 
Det a — c? and Dct/3 = \? are 5p(2,K) © Sy(2X) in- 
variant. Notice that, since either S'y(cr) or the quantities 
'S'y(c'"i)'s are invariant under local unitary operations, one 
has /(<t) = I {(Tsf)- Partial tracing of (Tsf over subsystem 
i yields <Ti = diag(a, a) and (T2 = diag(&, b), so that, 
finally, Eq. (I14> and Proposition 1 lead to Eq. (I26> . □ 
Notice that a and b constitute the symplectic spectra of, re- 
spectively, (Ti and (T2- Eq. (I26> emphasizes the relevant 
role played by the symplectic eigenvalues nzf{(Tsf) in de- 
termining the total amount of correlations contained in a 
quantum state of a continuous variable system, in striking 
analogy to the role played by the symplectic eigenvalues of 
the partial transpose of (Tsf in characterizing the amount of 
quantum correlations ll9l llOI . 

To better clarify this point, let us consider a symmetric 
state (Tsym, i-e. a state whose standard form fulfills a ^ b, 
so that its mutual information reads, according to Eq. i26\ 

licTsyrn) = 2f{a) - f [n ^ {(T syrn)] - / ["-+ (cTsym )] , 

(27) 

with symplectic eigenvalues = y/{a^ ci)(a =F C2). 
On the other hand, the symplectic eigenvalues of the 
partially transposed covariance matrix d-sym (obtained 
from (Tsym by switching the sign of C2, see |6|) are 
= n^{a-syrn) = ^/(a =F Ci)(a ± C2). In particular, 
for an entangled state, the smallest eigenvalue is 7t,_ = 
V(a-|ci|)(a-|c2|). ' 

The symplectic eigenvalue n_ encodes all the information 
about the entanglement of the state, since the necessary and 
sufficient criterion for entanglement reduces to n_ < 1/2, 
while the entanglement of formation Ep{asym) reads ^ 

EF{(Tsyni) = max{0,.g[n_]} , (28) 

iij|i,„(ii££)_ii#,.(it|a;) 

(see ||9l for details); it correctly reduces to /(cr)/2 for pure 

'if (T is entangled, then Det < 0, see Ref. fQ. 

^Ep{q) = rmn{p^,|y,.)} J] piE{\tpi)(tl)i\), v/here E{\^p){tp\) is the 
entropy of entanglement of the pure state | ■(/>), defined as the Von Neumann 
entropy of its reduced density matrix, and the min is taken over all the 
pure states realization of g = ^ pi \4>i)(i)i\- 



symmetric states. Even the quantification of quantum cor- 
relation provided by negativity ^8^, which is computable 
also for non symmetric states, reduces for a two-mode 
Gaussian state cr to a simple function of n_(cr). 

The dependence of entanglement on the eigenvalue 
n_(cr) is due to the fact that the biggest eigenvalue of 
the partially transposed covariance matrix a can be eas- 
ily shown to fulfill Ineq. i22\ . Thus alone can be 
responsible of the violation of the PPT ('positivity of the 
partial transpose') criterion for separability |6|, which can 
be recast as 

n-{(T) > i . (29) 

All the quantification of entanglement for two-mode Gaus- 
sian states available at present just quantify the violation of 
this inequality. 

We have extensively shown how both quantum and clas- 
sical correlations of a Gaussian state are encoded in sym- 
plectic spectra of the global, reduced and partially trans- 
posed covariance matrices of the state. 

5 Conclusions 

In conclusion, we have characterized mixedness and total 
correlations of two-mode Gaussian states by deriving their 
Von Neumann entropy and mutual information. Comparing 
these quantities with the entanglement of formation of sym- 
metric states shows that a crucial information about quan- 
tum and classical correlations lies in the symplectic eigen- 
values of the covariance matrix and of its partial transpose. 
The problem is still left open of determining a fully satis- 
factory quantification of the purely quantum correlations in 
a general two-mode Gaussian state. 
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